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, Abstract. Source blending in microlensing experiments is known to modify the 

Einstein time of the observed events. In this paper, we have conducted Monte- 



o 



X 



Carlo calculations, using the analytical relationships derived by Han (1999) to 
quantify the effect of blending on the observed event time distribution and op- 



, tical depth. We show that short-time events are affected significantly by source 

blending and that, for moderately blended sources, the optical depth r is globally 



overestimated, because of an underestimation of the exposure. For high blending 
situations, on the opposite, blending leads to an underestimation of the optical 
depth. Our results are in agreement with the most recent optical depth deter- 
minations toward the Galactic Center of the MACHO collaboration (Popowski 
et al. 2004) and the OGLE-H collaboration (Sumi et al. 2005) that use clump 
giants (less affected by the blending effect) as sources. The blending-corrected, 
lower optical depth toward the Galactic Bulge is now in good agreement with 
the value inferred from galactic models, reconciling theoretical and observational 
determinations. 



1. Introduction 

Microlensing is one of the most efficient ways to measure the amount of (sub)stellar 
matter, dark or luminous, in our Galaxy. By measuring the amplification of the flux of a 
source star, one can derive statistically the mass, distance and transverse velocity of the 
Send offprint requests to: Y. Alibert e-mail: yann. alibertSphim.unibe . ch 
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deflector. When a deflector passes near the line of sight between the observer and the 
source, the flux of the source is amplified by a factor (Paczyhski 1986) 



A{u) = — , „ where u{t) = -i/umin^ , , , 
uv w + 4 V V*E 

if the origin of time is taken at the moment of maximum amplification. The light curve 
depends on two parameters only, Umin and tE, where Umin is the minimum impact pa- 
rameter (the minimum distance to the line of sight, normalized to the Einstein radius 
of the lens Re), and tE is the Einstein time, i.e. the time required to cross the Einstein 
radius. 

The Einstein radius is given by the lens mass M, the observer-source distance L, and 
the ratio x of the lens-observer distance to L : 

-Re = — 5— ia;(l - x) 

When a microlensing event occurs, a theoretical light curve is ajusted to give Umin 
and Ie'- Wmin is related to the maximum amplification by Umin = /(^max), where the 
function / is defined by /(a;) = y^ ^Jf_^ — 2. The Einstein time is obtained from the 
event duration during which A > At (where At is the threshold amplification, usually 
taken as 2>l\f^, which corresponds to Ut = 1) by: 



„2 _ „ .2 



= id / \A(t" - Mmin 

The observed optical depth is estimated from the Einstein times of the events as 
(Alcock et al. 1997): 

Tobs = / ^(*Ei)) 

where E is the exposure, i.e. the number of observed sources times the duration of the 
experiment (in stars x years), and e is the detection efficiency. This optical depth can be 
compared with the one obtained from galactic models using 
47rG 



( 4ttG 
J ~ 



-x{l - x)L^p{xL)dx (1) 



for a constant observer-source distance, where p is the galactic mass density distribution. 

A long standing unsolved, major problem concerning the microlensing experiments 
towards the Galactic bulge is the significant discrepancy between the optical depths 
derived from microlensing experiments and the one calculated from theoretical models, 
as identified originally by Kiraga & Paczyhski (1994). Using Difference Image Analysis, 
and using main sequence stars as microlensing sources, the MACHO project (Alcock et 
al. 2000a) deduced an optical depth r_6 = 2.91^0 45 at 2a with 99 observed events, where 
T_6 = t/10~^. The OGLE collaboration derived t_6 = 3.3 ± 2.4 at 2a with 9 observed 
events (Udalski et al. 1994). More recently, the MOA collaboration has determined an 
optical depth t_6 = 2.59t[JJ^ (Sumi et al. 2003). 
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On the other hand, all Galactic models predict a typical optical depth r ~ 1 — 

1.2 X 10"^ (see for example Peale 1998, Bissantz et al. 1997). As shown by Peale (1998), 
one can reproduce the observed optical depth only if the bulge mass is equal to Mi, = 

3.3 X 1O^°M0, whereas gas dynamics calculations based on DIRBE observations to derive 
the gravitational potential of our Galaxy yield a maximum mass for the bulge Mt^^^ ~ 
2 X IO^^Mq (Bissantz et al. 1997, Englmaier et al. 1999, Sevenster et al. 1999). More 
recently, Han & Gould (2003) have derived an optical depth t_6 = 0.98, for events due 
to bulge stars, and T_g = 1.63 when including the events due to the disk. 

Using brighter stars (in particular red clump giants), the MACHO and the OGLE-H 
collaborations have derived recently new optical depths determinations, namely T_g = 
2.171^-3^ for 62 events for the MACHO group (Popowski et al. 2004), and r_6 = IMtoH 
for 33 events for the OGLE-II group (Sumi et al. 2005). These values, closer to the 
theoretical expectations, are significantly lower than the aforementioned ones based on 
main sequence star sources. This raises two important questions: 1) what is the origin of 
the systematic discrepancy between optical depth measured using faint (main sequence) 
and bright (clump giant) stars, and 2) which one should be compared with theoretical 
calculations. We demonstrate here that this difference is, for a large part, due to blending 
effects, and that the optical depth measured using main sequence stars is systematically 
biased. 

Microlensing experiments are subject to different biases like source blending (the 
observed lensed source is blended with other unresolved stars), lens blending (when the 
lens is a luminous star, its own luminosity may affect the observed light curve), and 
amplification bias (the amplified source is a star below the detection limit close to another 
star above the detection limit). Lens blending has been adressed by Han (1998); the effect 
is found to result in an underestimation of the real optical depth by ~ 10%. Amplification 
bias has been studied also by Han (1997): its effect is an overestimation of the optical 
depth by a factor ~ 1.71q'4, depending on the size of the seeing disk. Alard (1997) studied 
source blending by taking it into account directly in Eq.^ The effect on the optical depth 
is an overestimation of ^ 15%, but the effect on tE-histograms was not quantified. In this 
paper, we study in detail the afore-mentioned first kind of blending, i.e. the amplified 
source is a star above the detection limit, blended with a star above the detection limit. 

In the case of source blending, an observed source is the superposition of different 
source stars, all above the detection limit. When an event occurs, only one of the stars 
composing the observed source will be amplified, in general ^. The observed flux is no 
longer the baseline flux multiplied by the amplification, but instead, if one notes Fs the 
flux of the amplified source and Fq the fiux of the other unresolved sources, 

_ A{u{t))Fs + Fq 
^ the probability that more than one source is amplified is negligible. 
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It is very difficult to account correctly for this effect for, as shown by Han (1999), there 
is an almost perfect degeneracy between unblended and blended light curves. Moreover, 
it is very difficult to determine the amount of blending by photometrical means (see 
Wozniak & Paczyhski 1997). Nevertheless, using astrometric shift during a microlensing 
event, Goldberg and Wozniak (1998) have shown that about half of the sources observed 
by the OGLE collaboration suffer some amoimt of blending. 

To quantify the effect of source blending on the duration of a microlensing event and 
on the inferred optical depth, we first outline in §2 the relations between and Umin 
with and without blending. In §3, we derive the distribution hmction to be used for 
blending. The consequences on the Einstein time histograms and on the optical depth 
are examined in §4 and 5, respectively. Section 6 is devoted to the conclusions. 

2. Analytical relationships 

We will follow the previous calculations by Han (1999) with slightly different notations. 

When a star suffers blending, the total observed flux is the sum of the flux from the 
source Fs and the flux from the other unresolved stars Fq. We define the dimensionless 
parameters B as B = Fs/{Fs + Fq) and B = 1 — B. 

The observed amplification of a microlensing event thus reads: 

A{t) = B + ^ + B 

Using the observed light curve, one can derive some quantities, such as the impact pa- 
rameter, the duration of the event. We will label 'obs' the quantities derived when one 
supposes that B = 1 (i.e. when one ignores that the source is blended). The quantities 
corrected from blending will be called 'true', and will have no index. 
The observed impact parameter Umin,obs is 

^min,obs — f (^max,obs) 

= / (b "™^+ |^ + B 
SO that the duration of the event 



B 

is underestimated. 

We can then relate the observed Einstein time to the true Einstein time, the true 
impact parameter and the amount of blending by : 



*E,obs = teW P ^ ) ~ / \/ '"T^ - Mmin.obs^ (2) 



At- B 
B 

These relations allow to relate the 'observed' and 'real' quantities, for a known amount 
of blending. As shown by Han (1999), the curves characterized by (^e, Wmin, B) and the 
corresponding (iE.obs, Wmin.obs, -S = 1) lie very close to each other, and the difference is 
smaller than the observational errors. 
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3. Blending Distribution Function 

In order to estimate the effect of blending, one needs to know the distribution function 
(DF) of the blending factor B. We first need to estimate the number of stars whose pro- 
jection lies inside a small surface AS*. For that purpose, we suppose that the projected 
density of stars is constant and that the stars are distributed randomly. By doing so, we 
neglect all systematic effects in the spatial distribution of stars. Note that this approxi- 
mation is used by the MACHO group to compute their detection efficiency toward the 
LMC (Alcock ct al. 2001). Then, the probability that there are n stars inside AS is given 
by the Poisson law: p{n) = ri"e^"/n!, where n is the mean number of stars lying inside 
AS-. 

The second step is to calculate the probability that there are n stars inside the seeing 
disk, of surface AS", of an observed source. This probability P{n) is the probability that 
there are n stars inside AS, knowing that there is at least one star inside AS. We thus 
have 

Pin) ^ pin\n > 1) = = , f 

}_^,^j^p[i) n!(l-e ") 

The only free parameter in this law ^ is n, the mean number of stars inside a seeing disk. 

This parameter depends on the size of the seeing disk, and we will usually take n ~ 1.257, 

which ensures that half of the observed sources are unblended (i.e. P{1) = 0.5), as inferred 

by Goldberg and Wozniak (1998) from the OGLE data. Once the number n is given, Fq 

and Fs are randomly determined for a given bulge luminosity function (LF). 

For this latter, we use the LF obtained towards the Baade's Window with the Hubble 
Space Telescope (HST) (Holtzman et al. 1998). Since there are very few observed bright 
stars, we have proceeded as follows to obtain a combined LF: for bright stars {V < 18.5) 
the HST LF is matched to a power law function, with a —2 exponent. We restrict our 
calculations to sources fainter than Vsup — 16, which corresponds to the majority of 
observed sources (Alcock et al. 1997), and we performed calculations with different values 
of this parameter {Vsup ranging from 15.5 to 16.5). 

The faint cutoff Vinf is more difhcult to determine. We tried different limit magnitudes 
between V — 21.5 and V = 22.5, which correspond to the lower cutoff of the MACHO 
LF (Alcock et al. 1997). Since we are only concerned by source blending and not by 
amplification bias, there is no effect arising from the part of the LF below the detection 
limit. 



4. Effect on tE-histograms 

Blending affects the Einstein time histogram in two different ways. The first one is to 
decrease to tE,obs (Eq[21). The second one is to reduce the number of observed events: if 
^ In the following, we will refer to P as the Poisson law, although this is not exactly the 
standard Poisson law because of the different normalization factor 
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the maximum amplification is lower than the threshold amplification, the event is missed. 
This occurs if B is lower than -Bmin, where -Bmin is a function of Umin: 



The Einstein time histograms are calculated using a Monte-Carlo code originally 
developed by Chabrier and Mera (Mera, Chabrier & Schaeffer 1998). We have simulated 
10^ events which result, once including detection efhciency, in about ^ 10^ observed 
events. The Galactic model (model 1) used in the calculations is the following: a bar- 
shaped bulge with a density given by the model of Zhao (1996). The bar is in the plane 
of the disk, with an angle (j) — 13° with respect to the Sun-Galactic center line, the near 
end of the bar lying in the first quadrant. The bulge mass is taken to be 2.2 x 1O^°M0. 

The bar rotates rigidly with an angular speed Qha.r — 63.6km/s/kpc, the velocity 
dispersion is CTbuigc = 110 km/s in all directions. The model of the disk is a classical double 
exponential disk (Bahcall & Soneira 1980), with a scale length Rd — 2700 pc and a scale 
height h — 300 pc. The local normalization of the disk is pq = 0.05 Mq/pc^. The rotation 
velocity is Wrot,disk = 210km/s and the isotropic velocity dispersion is adisk = 20 km/s. 
The velocity distribution of each star is assumed to be gaussian, with the mean velocity 
equal to the rotation velocity and the afore-mentioned dispersions. The Sun is in the 
galactic plane, at a distance of Rq = 8 kpc from the Galactic center. Its velocity is equal 
to the one of the local standard of rest, Vq = ^lsr ~ 210 km/s. The visibility function 
of the sources is the one used by Kiraga & Pacsyhski (1994), with f3 — —1. The mass 
function is the one derived by Mera et al. (1998). 

We stress that the aim of the present paper is not to examine the validity of different 
galactic models, but to study the effect of source blending on the microlensing event time 
distribution and optical depth. 

To compute the time-histograms, we also need to take into account the detection 
efficiency. We used the clump giant efficiency toward the bulge: this efficiency is equal to 
the sampling efficiency and is not corrected for blending (Alcock et al. 1997). 

Figure n shows the effect of blending on the theoretical histogram (fraction of events 
as a function of the Einstein time) for different values of n in the Poisson distribution 
function of B (§3). As expected, the effect of blending is to decrease the mean Einstein 
time when the fraction of unblended stars / decreases. At the maximum of the histogram 
(between 4 and 6 days), there are 25% more events in the case of the Poisson DF {n ~ 
1.25) compared to the unblended histogram. This difference increases in the short-time 
region, and reaches 50% for — 2.5 days, but due to the decrease of the efficiency, 
there are very few observed events in this region. The mean Einstein time is 10% to 
15% lower for the three Poisson DF. If we convert this difference into a mean lens mass 
((<e) oc {y/m)), we obtain a difference larger than 20% between the unblended case and 
the Poisson law, i.e. the observed mean mass (derived from observations assuming that 
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Fig. 1. Einstein time histogram for differents amounts of blending. The curves give the 
relative number of short time events (below 66 days), binned in one-day time intervals. 
The curves are computed assuming a Poisson law with different values of n corresponding 
to different fractions of unblended stars (/). 

there is no blending) and the real mean mass (that coTild be derived from observations 
if one knows the amount of blending) i.e. (m)^^^ / (to)^^^! ^ 



5. The observed optical depth 

Source blending has three effects on the estimation of the optical depth r. The first one is 
to lower tg, which tends to underestimate the optical depth. The second one, as we just 
showed, is to underestimate the exposure (when a single source is thought to be observed, 
there are in fact many sources that may be lensed), which tends to overestimate the 
optical depth. The last one is due to the efficiency e which, depending on the variations 
of the efficiency with tB, can either underestimate or overestimate r. In order to quantify 
the resulting global effect, we performed the same Monte-Carlo calculations as described 
previously. 

We first need to relate the true exposure to the observed one. If n{i) is the real number 
of sources inside the observed source i, the ratio of the true to the observed exposure is 
given by: 

E ^ Esourccs"(0 

^ohs ^ sources 

where A^sources is the number of monitored sources. The true exposure is given by 

oo 

E=Y,nx N{n) 
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where N{n) = -/Vgourccs x ^(") is the number of observed sources composed of n blended 
sources. Therefore 

n>l 

It is then straightforward to compute the ratio of the observed optical depth to the 
real optical depth, for a given DF of the blending B. The luminosity of the source is 
calculated for each event with the bulge LF. Then, knowing the real number n(i) of 
sources, we can calculate the amount of blending. The ratio of the true to observed 
optical depth thus reads: 

T _ £-obs \ ^ / \ ^ ^E,obs 

Tohs ^ E ^ £(<e) / ^ e(iE.obs) 

Figure [3 shows the ratio r/robs for different values of the fraction of unblended stars 
/. / is related to the parameter of the Poisson law by 

1-e-" ^ ' 

In order to test the sensitivity of this ratio to the Galactic model, we have performed the 
same calculations for two other models. Model 2 is obtained by changing only the bar 
angle (set to (f) — 20°), the mass of the bulge (1.2 x \Q^^Mq), and the disk scale length 
[Rd = 3500 pc), while the local normalization of the disk pq is unchanged. Model 3 is 
obtained by changing the velocity parameters: (iTbuigo, Cdisk) — (100, 30) km/s, ilbar = 
50.0km/s/kpc and Urot.disk = 180 km/s. As seen on Figure 2, the ratio is not very 
sensitive to the Galactic model. This was expected since, forgeting the e term, tE,obs = 
iE X g{unun,B) with 



^min 



g (Wmin, B) = , „ ^ 6 {B - Smin) 

^ UT - timin,obs 

where is the Heaviside function (equal to for a negative argument). Then the ratio 
Tobs/T can be written: 

Tobs/r — - {g (Umin, B)) , 

1 — e " 

which is independent of the galactic model. The ratio is also quite insensitive to the cutoff 
values of the LF since all the curves in the case of model 1 are very close to each other 
and the differences are of the order of the errors in the Monte-Carlo calculations. 

The effect of blending is to overestimate the optical depth, for each of the three models 
considered. The ratio of the true optical depth to the observed one can reach 75%. This 
result is consistent with the one derived by Alard (1997) using a different mehod. By 
taking into account the blending directly in Eq. ^ he found an overestimation of ~ 15%, 
whereas we obtain ^ 20% for / = 0.5. 

We can relate our results to the last determinations of the optical depth using 
clump giant stars as sources. Popowski et al. (2004) argue that, using clump giant and 
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Fig. 2. Ratio r/robs as a function of the unblended source fraction for different values of 
the LF cutoff. Standard values are \^nf = 22 and V^up = 16 and these values are assumed 
except when indicated by the legend. Triangles and circles are computed using the same 
LF as the continuous line, but with model 2 and model 3 (see text). 

Differential Image Analysis, they can derive an optical depth with 110 blending effect. 
This comes from the fact that, dne to their luminosity, clump giant stars are unlikely 
to suffer large amounts of blending. Note, however, that Sumi et al. (2005) suggest that 
even clump giant stars can suffer some amount of blending, at least in the OGLE-II 
sample. Using the last MACHO value as the true optical depth, the corresponding ob- 
served optical depth for blended sources, assuming a Poisson law with / ~ 0.5, would 
be T_6 = 2.71^0 47, compatible with the value obtained from main sequence sources, 
T_6 = 2.9i;^:4^ (Alcock et al. 2000). 

Note that for very crowded fields, n > 10, the ratio Tobs/''' is greater than one, so that 
in that case blending effects yield an underestimation of the optical depth. 

6. Conclusion 

We have shown in this paper that blending can have an important effect both on the 
Einstein time histograms and the inferred optical depth. 

For histograms, the effect is to lower the observed t^ and then to decrease the 
derived mean lens mass. The change in the histogram is about 15% in the low-tE region 
corresponding to the maximum of the expected histogram, and can decrease substantially 
the inferred mean lens mass (the difference can reach 20%). 
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For the optical depth, wc have shown that, for all values considered for n, the effect 
is to overestimate the optical depth as a result of an underestimation of the exposure. In 
the case of the Poisson law we used, with an unblended fraction aroimd 50%, this effect 
can reach 20%. These results arc quite insensitive to the cutoff of the LF and to the 
galactic model. We note, however, that in the case of very high blending, for n > 10, the 
effect of blending is to underestimate the optical depth. This can occur for microlensing 
experiments towards very crowded fields, like M31 (AGAPE project, Ansari et al. 1997 
and references therein). 

Finally, to relate the observed optical depth to the one derived from different models, 
it is necessary to combine these results with the effect of other biases, like lens blending 
(Han 1998) and amplification bias (Han 1997, Alard 1997). The effect of lens blending 
is to underestimate the optical depth by about 10%, which can cancel, at least in part, 
the effect of source blending derived presently in the case of moderate blending (for the 
Poisson law used in these calculations, this corresponds to an unblended fraction larger 
than 80%). For amplification bias, Han (1997) has shown that the effect would be to 
overestimate the optical depth. Therefore, the net effect of the various blending effects 
is very likely to overestimate the correct optical depth by a substantial fraction. Taking 
blending effects into account brings into agreement the optical depth derived with main 
sequence sources (Alcock et al. 2000, Udalski et al. 1994), the one derived using brighter 
stars (Popowski et al. 2004, Sumi et al. 2005), and the one derived from usual Galactic 
models (Bissantz et al. 1997, Englmaier et al. 1999, Sevenster et al. 1999) at the 2cr 
level, and thus reconciles experimental and theoretical determinations without drastic 
changes in galactic modelling. Further observational determinations of the exact amount 
of blending are needed to nail down precisely the net effect of blending. 
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